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ABSTRACT 

We present the derivation of two-integral distribution functions for the first four mem- 
bers o f the family of generalized Kalnajs discs, recently obtained by I Gonzalez fc Reinal 
( 2006), and which represent a family of axially symmetric galaxy models with finite 
radius and well behaved surface mass density. In order to do this we employ several 
approaches that have been developed sta rting from the p otential-density pair and, es- 
sentially using the method introduced bv lKalnali (|1976f) . we obtain some distribution 
functions that depend on the Jacobi's integral. Now, as this method demands that the 
mass density can be properly expressed as a function of the gravitational potential, 
we can do this only for the first four discs of the family. We also find another kind 
of distribution functions by starting with the even part of the previous distribution 
functions and using the maximum entropy principle in o rder to find the o dd part and 
so a new distribution function, as it was pointed out bv lDeion ghc (198(|). The result 
is a wide variety of equilibrium states corresponding to several self-consistent finite 
flat galaxy models. 

Key words: stellar dynamics - galaxies: kinematics and dynamics. 



1 INTRODUCTION 

The problem of finding self-consistent stellar models for 
galaxies is of wide interest in astrophysics. Usually, once the 
potential-density pair (PDP) is formulated as a model for a 
galaxy, the next step is to find the corresponding distribu- 
tion function (DF). This is one of the fundamental quantities 
in galactic dynamics specifying the distribution of the stars 
in the phase-space of positions and velocities. Although the 
DF can generally not be measured directly, there are some 
observationally accesible quantities that are closed related 
to the DF: the projected density and the light-of-sight ve- 
locity, provided by photometric and kinematic observations, 
are examples of DF's moments. Thus, the formulation of 
a PDP with its corresponding equilibrium DFs establish a 
self-consistent stellar model that can be corroborated by as- 
tronomical observations. 

On the other hand, a fact that is usually assumed in as- 
trophysics is that the main part of the mass of a typical spi- 
ral gal axy is concentrated in a thin disc (jBinnev fc Tremainel 
(|l987T ) V Accordingly, the study of the gravitational poten- 
tial generated by an idealized thin disc is a problem of great 
astrophysical relevance and so, through the years, different 
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approaches has been us ed to obtain th e PDP for suc h kind 
of th in disc models (see iKuzm in (1956) and Toomre (|l963l : 
1 1964T ) . as examples). Now, a simple method to obt ain the 
PDP of thin discs of finite radius was developed by iHunterl 
the simplest exa mple of d i sc obt ained by this method 
being the Kalnajs disc (|Kalnaisi (|l972t 0. 

In a previous paper (|Gonzalez fc Reinal (2006)) we 
use the Hunter method in order to obtain an infinite 
family of axially symmetric finite thin discs, character- 
ized by a well-behaved surface density, whose fir st mem- 
ber i s precisely the well-known Kalnajs disc l|Kalnaisl 
fl972)) and, more recently, a detailed study of the 
kinematics characterizing such discs was be made by 
iRamos-Caro. Lopez-Suspes fc Gonzalez! (|200&f ) . So, in order 
to obtain self-consistent stellar models for this family of gen- 
eralized Kalnajs discs, we will consider at the present paper 
the derivation of some two-integral DFs for the first four 
members of the family. 

Now, as is stated by Jeans's theorem, an equilibrium 
DF is a function of the isolating integrals of motion that 
are conserved in each orbit and, as it has been shown, it 
is possible to find such kind of DFs for PDPs such that 
there is a certain relationship between the mass density and 
the gravitational potential. The simplest case of physical 
interest corresponds to spherically symmetric PDPs, which 
are described by isotropic DFs that de pends on the tota l 
energy E. Indeed, as was be shown by lEddingtonl l|l916h . 
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it is possible to obtain this kind of isotropic DFs by first 
expressing the density as a function of the potential and 
then solving a n Abel integral equation. 

Recently, Ijiang fc Ossipkovl (|2007l ) found that a simi- 
lar procedure can be performed in the axially symmetric 
case, where the equilibrium DF depends on the energy E 
and the angular momentum about the axis of symmetry L z , 
i.e. the two classical integrals of motion. They developed 
a formalism that essentially co mbines both the Eddington 
formulae and the IFrickd l| 19521 ) expansion in order to ob- 
tain the DF's even part, starting from a density that can 
be expressed as a function of the radial coordinate and the 
gravitational potential. Once such even part is determined, 
the DF's odd part can be obtained by introducing some 
reasonable assumptions about the mean circular velocity or 
using the maximum entropy principle. 

On the other hand, another method appropriated to 
find DFs depending only of the Jacobi's integral, E r = 
E — QL Z , fo r axially symm etric flat galaxy models was in- 
troduced bv lKalnaia (|l976l ') for the case of disc- like systems 
and basically consists in to express the DF as a derivative 
of the surface mass density with respect to the gravitational 
potential. Such method does not demands solving an inte- 
gral equation, but instead is necessary properly to express 
the mass density as a function of the gravitational potential, 
a procedure that is only possible in some cases. 

In this paper we will use both of the mentioned ap- 
proaches in order to obtain equilibrium DFs for some gener- 
alized Kalnajs discs. Accordingly, the paper is organized as 
follows. First, in Section [2] we present the fundamental as- 
pects of the two methods that we will use in order to obtain 
the equilibrium DFs. Then, in Section [3] we present a sum- 
mary of the main aspects of the generalized Kalnajs discs, 
at Section 13.11 and then, in Sections 13.21 - 13.51 we derive 
the DFs for the first four members of the family. Finally, in 
Section 21 we summarize our main results. 



2 FORMULATION OF THE METHODS 

Assume that $ and E are, respectively, the gravitational 
potential and the energy of a star in a stellar system. One 
can choose a constant <3>o such that the system has only 
stars of the energy E < $0, and then define a relative po- 
tential = — $ + <& q and a relative energy e = —E + $0 
l|Binnev fc Tremaind (|l987l )). such that e = is the energy 
of escape from the system. Both the mass density p(r) and 
the DF /(r, v) (r is a position vector and v is a velocity 
vector) are related to ^(r) through Poisson's equation 



V 2 * = AixGp = AtiG I /d 3 v, 



(1) 



where G is the gravitational constant. 

For the case of an axially symmetric system, it is cus- 
tomary to use cylindrical polar coordinates (R, <p, z), where 
v is denoted by v = (vr, %>,w z ). As it is well known, 
such system admits two isolating integrals for any orbit: 
the component of the angular momentum about the 2-axis, 
L z = Rv v , and the relative energy e. Hence, by the Jeans 
theorem, the DF of a steady-state stellar system in an axi- 
ally symmetric potential can be expressed as a non-negative 
function of e and L z , denoted by f(e,L z ). Such DF, that 



vanishes for e < 0, is related to the mass density through 
eq. ©. 

In this subject, f(e,L z ) is usually separated into even 
and odd parts, f+(e,L z ) and f-(e,L z ) respectively, with 
respect to the angular momentum L z , where 

/+(£,£,) = §[/(e,L.) +f(e,-L,)], (2) 

and 

f-{e,L z ) = \[f{e,L z )-f(e,-L z )]. (3) 

So, by using s = <P — \ ( v% + vf„ + v%), the integral given by 
{TJ can be expressed as l|Binnev fc Tremaind (119871 )) 



47T 

P= R 



f+{e,L z )dL z 



ds. 



(4) 



For a given mass density, this relation can be considered as 
the integral equation determining f+(e,L z ), while the odd 
part satisfies the relation 



pR(v v ) 



47T 



-fV 2 (*- £ ) 



L z f-(e,L z )dL z 



de. (5) 



This integral equatio n was fi rst fou nd bv lLvnden-Belil (|l962h 
and then applied by lEvansl fl993) into calculating the odd 
DF for the Binney model, under the assumption of {v v ) 
having some realistic rotational laws. 

As is not known, we cannot compute f-(e,L z ) 

directly by eq. (JSJ but what we can do is to obtain the 
most probable distribution functions under some suitable 
assumptions. Once f+(e, L z ) is known, /_ (e, L z ), and there- 
fore /(e,L 2 ), can be obtained b y me ans of the maximum 
entropy principle (see iDeionghd 1 19861) ) and we obtain 

„ aL ' — 1 



f-(e,L z ) =/ + ( £ , L z )- 



f(s,L z ) = 



2f+(e,L z ) 
1 + e - aL * ' 



(6) 



(7) 



where a is the parameter depending on the total angular 
momentum. Obviously, \f-(e, L z )\ < f+(e,L z ). Also, the 
system is non-rotating when a — and maximally rotat- 
ing as a — + 00, i.e., for a — > +00, it is anticlockwise and 
f(s y L z ) — [1 + sign(L z )]/+(e, L z ), for a — » —00, the rota- 
tion is clockwise and f(e,L z ) = [1 — sign(L z )]/+(e, L z ). The 
parameter a reflects the rotational characteristics of the sys- 
tem. 

As it was point e d out by IFrickd l|l952l ) and recently 
bv lJiang fc Ossipkovl l|2007h . the implementation of integral 
equations @ and demands that one can express p as a 
function of R and This holds indeed for the case of disc- 
like systems, which surface mass density E is related to / 
through 



E = 4 



f + (s y L z )dL z 



V2i? 2 (*-e)-Ll 



de. 



(8) 



In order to incorporate the formalisms developed for the 
3-dimensional case to deal with disc-like systems, we have 
to gen erate a pseudo-volume density p l|Hunter fc Quiar] 
(1993)), according to 



Jo V*- *' : 



(9) 
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which must take the place of p in Q and ([5]). In particular, 

m 

when p(^!,R) — Yl R 2n Qn{^), the corresponding even DF 
is (jjiang fc Ossipkovl (|2007T O 



f+{e,L z ) 



2 — n — 3/2j^2n. 

~ J) 



7r 3/2 r ( 



J_ /cT+V(£)\ 

V rf* n+1 A 



(10) 



Another simpler method to find equilibrium DFs cor- 
responding to axially symm etric disc- like systems, was in- 
troduced bv lKalnaisi (|l976l l. Such formalism deal with DFs 
that depend on the Jacobi's integral E T = E — Q.L Z , i.e. the 
energy measured in a frame rotating with constant angular 
velocity Q (subscript r denotes a quantity measured in the 
rotating frame). It is convenient to define an effective po- 
tential = <J? — ^Q 2 R 2 in such way that, if we choose a 
frame in which the velocity distribution is isotropic, the DF 
will be L z -independent and, from ([5J, the relation between 
the surface mass density and the DF is reduced to 



E = 2tt 



f{e r )de r - 



(11) 



Here, * r = \Q?R 2 + $ r and e r = e + QL z +$ r - <&o, 
i.e. the relative potential and the relative energy measured 
in the rotating frame. Moreover, if one can express E as a 
function of differentiating both sides of (|ll|l with respect 
to <]/,., we obtain 



/(£r) 



1 dE 

2tt d* r 



(12) 



Note that in this formalism it is also necessary to express 
the mass density as a function of the relative potential. 



3 DFS FOR THE GENERALIZED KALNAJS 
DISCS 

3.1 The generalized Kalnajs discs 

The generalized Kalnajs discs fo rm an infinite family of ax i- 
ally symmetric finite thin discs (jGonzalez fc Reinal ([2006)). 
The mass surface density of each model (labeled with the 
positive integer m > 1) is given by 

?2 -l m-1/2 



E m (i?) = E 



with the constants Ec given by 

{m) _ (2m + l)M 



(m) 



1-^ 



(13) 



(14) 



2na 2 ' 

where M is the total mass and a is the disc radius. Such 
mass distribution generates an axially symmetric gravita- 
tional potential, that can be written as 



3\n (£,?)) =-J2C 2 nq2 n (0 P 2n(v) 



(15) 



where P2n(r]) and <72n(£) = i 2n+1 Q2n(i£,) are the usual Leg- 
endre polynomials and the Legendre functions of the second 
kind respectively, and Gn n are constants given by 

MG*7r 1/2 (4n + l)(2m + 1)! 



C* 2 



a2 2m + 1 (2n + l)(m - n)W(m + n + |)g 2 n+i(0) ' 



where G is the gravitational constant. Here, — 1 < 77 < 1 and 
< ^ < 00 are spheroidal oblate coordinates, related to the 
usual cylindrical coordinates (R, z) through the relations 

R 2 = a + rf), z = ar)Z. (16) 

In particular, we are interested in the gravitational potential 
at the disc, where z = and < R < a, so£ = and 
77= ^l-R 2 /a 2 . 

If we choose "I>om in such a way that ^ m > 0, the cor- 
responding relative potential for the first four members will 
be 

^{R) = ^-{a 2 -R% (17a) 

*2{R)= l ^^^-m 2 a 2 + m\ (i7b) 

v &3 ( jR) = 5g^(ii a 6 _ 2AR 2 a 4 + 18R 4 a 2 - 5R 6 ), (17c) 
MR) - 5S??(93« 8 ~ 256i? 2 a 6 + 288i? 4 a 4 
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160R 6 a 2 + 35R 8 ), (17d) 



We shall restrict our attention to these four members. The 
formulae showed above defines the relative potentials that 
will be used to calculate the DFs by the implementation of 
the methods sketched in section [2] 



3.2 DFs for the m = 1 disc 

Given the mass surface distribution (|13p and the relative 
potential (|17a|) . we can easily obtain the following relation: 



Ei(*i 



±/2S« /2 

r\ 1 ' 

Had 



(18) 



where n = [37rGM/(4a 3 )] 1/2 . As we are dealing with disc- 
like systems, it is necessary to compute the pseudo- volume 
density p by eq. ([9]) in order to perform the integral Q, 



7rE (1> 

pi(*i) = -Fnr*i- 

1 LqO, 



(19) 



By using the right part of eq. (|10|l with n = 0, which is equiv- 
alent to calculate the Fricke component of (|19[). we obtain 
the even part of a DF which depends only on the relative 
energy 



fiiHe) 



27rr2oav2£ 



(20) 



At this point, we may no tice that this fi+(e) corre sponds 
to the DF formulated by iBinnev fc Tremaind (|l987i ) when 
fl = 0. (Note that there is a difference of constants as a 
result of a different definition of the relative potential and 
relative energy.) 

To obtain a full DF, we use the maximum entropy prin- 
ciple by means of eq. J7J and the result is 



ft A) (e,L z 



— • (21) 

As it is shown in figure [1] a determines a particular ro- 
tational state in the stellar system (from here on we set 
G — a = M = 1 in order to generate the graphics, without 
loss of generality). As a increases, the probability to find 
a star with positive L z increases as well. A similar result 
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Figure 1. Contours of f\ A \ given by l|21j l. for (a) a = 10; (b) 
a = 1; (c) a = 0.1; (d) a = 0. Larger values of the DF corresponds 
to lighter zones. 



Figure 2. Contours of f\ B \ given by l|26|l . for (a) Q = — 7r/4; 
(b) f2 = — 7r/16; (c) Q = 7r/4; (d) H = 7r/16. Larger values of the 
DF corresponds to lighter zones. 



can be obtained for a < 0, when the probability to find 
a star with negative L z decreases as a decreases, and the 
corresponding plots would be analogous to figure [1] after a 
reflection about L z — 0. 

We can generalize this result if we perform the analysis 
in a rotating frame. At first instance, it is necessary to deal 
with the effective potential in order to take into account 
the fictitious forces. Choosing conveniently $or, the relative 
potential in the rotating frame takes the form 



(«0 



'(0 



R 2 ), 



(22) 



so the corresponding mass surface density and the pseudo- 
volume density can be expressed as 

/n V (l) 

si(*iQ= *t /2 (23) 



ly/Q,l - Q 2 



and 



p(tfir-) = 



ttE 



(24) 



ls/% - n 2 

The resulting even part of the DF in the rotating frame is 



si 1 ' 



(25) 



and it can be derived following the same procedure used to 
find (e) or by the direct application of eq. (|12p . 

Finally, one can come back to the original frame through 
the relation between e r , e and L z to obtain 



fl B \e,L z 



[2(e + flL z ) - Q 2 a 2 
2-nas/Q.l - Q 2 



-1/2 



(26) 



which is totally equivalent to the Binney and Tremaine DF 
for the Kalnajs disc. Its contours, showed in figure[2j reveals 
that the probability to find stars with e < Q 2 a 2 /2 — QL Z is 



zero, while it has a maximum when e > Q 2 a 2 /2 — QL Z (this 
defines the white strip shown in the figure) and decreases 
as e increases. Nevertheless, this f[ B \e,L z ) is quiet unre- 
alistic as it represent a state with {v v ) = QR and it means 
that the behavior of the system, in disagreement with the 
observations, behaves like a rigid solid. 

However, we can generate a better DF if we took only 
the even part of (I26|) and using the maximum entropy prin- 
ciple to obtain 



,L Z ) = 



Vil\e,L z 



(27) 



1 + e~ aL * 

As it is shown in figure [3j there is a zone of zero proba- 
bility, just in the intersection of black zones produced by 
f{ B (e,±L z ), and there are also two maximum probability 
stripes. The variation of the Q parameter leads to the change 
of inclination of the maximum probability stripes and it is 
easy to see that the DF would be invariant under the sign of 
Q, by the definition of the even part. Furthermore, a plays 
a similar role than in figure [T] increasing the probability of 
finding stars with high L z as a increases and vice versa. 

3.3 DFs for the m = 2 disc 

Working in a rotating frame we found that the relative po- 
tential is given by 



R 



128a 5 



[45GM7r(a 2 - R 2 ) 
+30GMna 2 - 64a 



(28) 



while the surface mass density is given by (|13|l when m = 2. 
This case is a little more complicated than the usual Kalnajs 
disc, because the analytical solution of the pseudo-volume 
density cannot be performed with total freedom. We will 
need to operate in a rotating frame conveniently chosen in 
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Figure 3. Contours of f^, given by l|27|l . for Q = tt/4 and (a) 
a = 1; (b) a = 10. For C = tt/16 we have (c) a = 1; (d) a = 10. 
Larger values of the DF corresponds to lighter zones. 



such a way that the relation between the mass surface den- 
sity and the relative potential becomes simpler. 

From eq. (|28[) it is possible to see that if we choose the 
angular velocity as 



n = ± 



15GMtt 
32a 3 ' 



the relative potential is reduced to 

*2r(i?) = 



45GM7T, 2 D 2x2 

(a — R ) . 



(29) 



(30) 



128a 5 

Now, we can express the surface mass density easily in terms 
of the relative potential by the relation 

3/4 

-*2r ) , (31) 



E 2 (* 2r ) = E< 2 > ( 



\45GMty 

and the integral for the pseudo-volume density can be per- 
formed and we obtain 



P(*2r) 



2 7 rEr , r(7/4) / 128a 



T(9/4) \45GM-k 



3/4 



5/4 



(32) 



Using then eq. (|10p . the resulting even part of the DF 
in the rotating frame and the full DF in the original frame 
are 



and 



f ( 2 A \e,L z ) = K (e + OL z - §fiV) 



2 2\-i/4 



(33) 



(34) 



respectively, where Q. could take the values according to (29 
and k is the constant given by 



3Ei 2) / 128a 



3/4 



(35) 



8tt \45GM-k / 
Moreover, using the same arguments given in section 1)3.30 . 



Figure 4. In (a) and (b) we show the contours of f 2 A ^ , given by 
(3D, for = ^15tt/32 and fi = --^15^/32, respectively. The 
contours of f 2 A \ given by (36}, are plotted in (c) a = 1 and (d) 
a = 10. Larger values of the DF corresponds to lighter zones. 



it's convenient to take the even part of ()34)l and obtain a 
new DF using eq. (01, which is given by 



;^( e ,L,)- 2/2 + )(e ' L * 



l + e- 



(36) 



A more general case can be derived without the assump- 
tion (|29|) if we use the Kalnajs method in order to avoid the 
pseudo-volume density integral. Here we work in terms of 
the spheroidal oblate coordinates to obtain more easier the 
relation between the mass surface density and the relative 
potential, which can be expressed as 



T , , 45GMtt 4 15GAIit — 32a 3 f2 2 2 

* 2 " (??) = ^28^ + 64a " 

Now, it is possible to rewrite this expression as 

*2r(j?) = (KIT? 2 + K2) 2 + K3, 

with 



«1 



45GMtt 
128a : 



K2 = 



128a 15GMtt - 32a 3 £T 



45GMtt 



128a 



and 



«3 = 



128a /15GM-K - 32a 3 n 



3r>2 \ 2 



(37) 
(38) 

(39) 
(40) 

(41) 



45GMtt V 128a 
Then, as E can be expressed in terms of rj in the form 
E m fa) = E<"V m -\ (42) 
the relation between E 2 and ^2r is 



6 J. F. Pedraza, J. Ramos- Cam and G. A. Gonzalez 




( B) 

Figure 5. In (a) and (b) we show the contours of /j given by 
II45I I. for Q = n/10 and Q = ir/5, respectively. The remaining 
figures exhibit the contours of ft B , given by ]46t . for f2 = w/10 
and (c) a = 0.1; (d) a = 5. Larger values of the DF corresponds 
to lighter zones. 



Now, by using eq. (|12[1 . we obtain 

Ve7 



k 3 - k 2 



-,1/2 



£r - k 3 

and the result in the original frame is 
3E C 



(44) 



87TfcJ 



3/2 



^e + Q.L z -Q. 2 a 2 /2-k 3 



e + QL z -Q. 2 a 2 /2- fa 



(45) 

Obviously, this DF is the same as ([34} when the condition 
(f2"9")l is satisfied. Finally, by eq. @, the resulting DF with 
maximum entropy is given by 



f(B) ( j- s _ 2/ 2 ( + ) (e,£ z ) 



(46) 



We can see the behavior of these DFs in figures [4] and [5] 
In figures[4][a) and[4jb) we show the contours of for the 
two rotational states given by (I29|l . Such DF is maximum 
over a narrow diagonal strep, near to the zero probability re- 
gion and, similarly to the case showed in figure [3] the proba- 
bility decreases as e increases. Stellar systems characterized 
by different f2 are shown in figures [5j a) and[5jb), where we 
plot the contours of (this equals to f% when Q is given 
by (|29[0 . In this case the DF varies more rapidly as fl de- 
creases, originating narrower bands. The remaining figures 
show the contours of ft an d ft f° r different values of 

~( B) 

the parameter a, showing a similar behavior than ft . 



3.4 DFs for the m = 3 disc 

Once again, if we want to use Kalnajs method, it is necessary 
to derive the relation £3(^3,.) and, according to (|42l) . it is 



posible if we can invert the equation of the relative potencial, 
which in this case is given by 



T . . 175GMtt 6 105GMtt 4 

*3r(??) =^T7T ?7 + — — Tj 



512a 



512a 



(47) 



+ 



512a 



in order to obtain r)(^ 3r ). To solve it, we must deal with a 
cubic equation and with its non-trivial solutions; fortunately, 
we still have il as a free parameter. One can easily note that 
it is possible to write (|47p as 



*3r(??) = {l^lll 2 + K 2) 3 + K 3 , 



with 



Ki 



K2 



/175GMtt 
V 512a 

_ / 7ttGM \ 
~ V 2560a ) 
7GMn 



1/3 



1/3 



K3 



2560a ' 
and Q has be chosen as 



Q = ± 



21GMtt 
64a 3 ' 



Now, by replacing (|42[) into (148 1) . we obtain 



(3) 



E 3 =S 



and, by using eq. (|12[1 . 



(#3e ~ K 3 ) 1/3 - K 2 
Kl 



5/2 



5E 



(3) 



/s(er) = 



((Er - ^3) V3 - fSa) 



3/2 



127r Kl /2 (e r -K 3 ) 2/3 
Coming back to the original frame, the result is 



h{e,L z ) = 



5Ei 3) ((e + Q,L Z - QV/2 - k 3 ) 1/3 - « 2 ) 



127TK! 72 (e + OL z - fi 2 a 2 /2 - k 3 ) 



2/3 



while the respective DF with maximum entropy is 

2/3+ (e,L.) 



1 + e- 



(48) 

(49) 

(50) 
(51) 

(52) 
(53) 

(54) 

3/2 

(55) 
(56) 



In figure [6f a) we show the contour of / 3 , while in figures 
Ob))O c ) and[5{d) are plotted the contours of ft for differ- 
ent values of a. We can see that the behavior of these DFs 
is opposite to the previous cases. As the Jacobi's integral 
icreases, the DF also increases. 



3.5 DFs for the m — 4 disc 



As we saw in section 13.41 in order to find a DF using the 
Kalnajs method, we must find rjasa function of the relative 
potential in order to obtain E(^ r r .). For the m — 4 disc, the 
relative potential can be expressed as 



11025GA/7T 8 1575GMtt 
*4r = „„„„„ rf + 



32768a 



8192a 



+ 



2835GMtt 4 1575GMtt - 4096a 3 Q. 2 



(57) 



16384a 



-rf + 



8192a 
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Figure 6. In (a) we show the contours of /3 given by l|55| l and 
taking the positive value of Q. The remaining figures exhibit the 
contours of fz, given by ( 1561 . for (b) a = 10; (c) a = 1; (d) 
a = 0.1. Larger values of the DF corresponds to lighter zones. 



Figure 7. In (a) we show the contours of /4 given by l|66| l and 
taking the positive value of Q. The remaining figures exhibit the 
contours of /j, given by H68I I. for (b) a = 10; (c) a = 1; (d) 
a = 0.1. Larger values of the DF corresponds to lighter zones. 



Now, although we have to deal with a quartic equation, it 
is possible to rewrite (|57[) as 



where 



where 
«i 

K2 



K3 = 



K'4 



*4r = [(KIT] 2 + K 2 ) 2 + + «4, 

/ 1L025ttGM\ 1/4 
V 32768a J 

/ 225ttGM \ 1/4 
V 1605632a / 

81ttGM \ 1/2 
6272a / 
25281?rGM 



1605632a ' 
and f2 must to be chosen as 

n = ± 



135GM7T 
448a 3 ' 



(58) 

(59) 

(60) 

(61) 
(62) 

(63) 



Finally, by using (|42|l and (148 1) . we find the expression 



4 r — K4 — K3 — K2 



7/2 



(64) 



which, by means of eq. I|12[l. can be used to derive the even 
part of the DF in the rotating frame, 



7E 



(4) 



/4(£r) = 



K-4 — K3 — K2 



5/2 



16ttKj / ' 2 \/ (V£r — K4 — K 3 )(e r — re 4 ) 



(65) 



Therefore, the corresponding DF in the original frame is 
given by 



fi{e,L z ) = 



7Ei 4) [ 5 (e,L z )- K2 



18/2 



16ttkI /2 'g(e, L z )[g 2 (e, L z ) + k 3 ] 



(66) 



g(e, L z ) = \J a/s + OL z - n 2 a 2 /2 - k 4 - « 3 , (67) 
and the respective DF with maximum entropy is given by 

2/4+ (e,£.) 



1 + e- 



(68) 



In figure [7{ a) we show the contour of /4, while in figures 
[7Jb),[7jc) and[7Jd) are plotted the contours of / 4 for different 
values of a. As we can see, the behavior is analogous to the 
showed at figure [S] 



4 CONCLUDING REMARKS 

We presented the derivation of two-integral equilibrium DFs 
for some members of t he family of general i zed K alnajs discs 
previously obtained bv lGonzalez &: Reinal l|2006h . Such two- 
integral DFs were obtained, esentially, by expresing them 
as functionals of the Jacobi's integral, as it was sketched 
in the formalism developed by iKalnaisI (|l976l ) . Now, since 
such formalism demands that the surface mass density can 
be written as a potential-dependent function, the above pro- 
cedure can only be implemented for the first four members 
of the family, the discs with m — 1, 2, 3, 4. 

Indeed, the procedure requires that the expression given 
the relative potential *& r as a function of the spheroidal vari- 
able 77 can be analytically inverted in order to express the 
surface mass density E as a function of the relative potential. 
So, we can do this in a simple way for the discs with m = 1 
to 4. However, when m > 4 we must to solve an equation 
of grade larger than four, whose analytical solution do not 
exists. 

For the first two members of the family, the discs with 
m = 1,2, we also use the method introduced recently by 
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Ijiang fc Ossipkovl (|2007l ) in order to find the even part of the 
DF and then, by introducing the maximum entropy princi- 
ple, we can determines the full DF. This procedure was also 
used for the other three discs, starting from the Kalnajs 
method, so defining another class of two-integral DFs. Such 
kind of DFs describes stellar systems with a preferred rota- 
tional state, characterized by the parameter a. 

This paper can be considered as a natural comple- 
ment of the work previously presented b y Gonzal ez fc Reina 
|2006t ) andU amos-Caro. Lopez-Suspes fc Gonzalez! (|2008h . 
where the PDP formulation and the kinematics, respectively, 
of the generalized Kalnajs discs were analyzed. Now, by the 
construction of the corresponding two-integral DFs, the first 
four members of this family can be considered as a set of 
self-consistent stellar models for axially symmetric galaxies. 
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